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Abstract 

,The impact of pairing correlations on the fission barriers is investigated in Relativistic Hartree Bogoliubov (RHB) theory and 
Relativistic Mean Field (RMF)+BCS calculations. It is concluded that the constant gap approximation in the usual RMF+BCS cal- 
culations does not provide an adequate description of the barriers. The RHB calculations show that there is a substantial diff'erence 
in the predicted barrier heights between zero-range and finite range pairing forces even in the case when the pairing strengths of 
these two forces are adjusted to the same value of the pairing gap at the ground state. 
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1. Introduction 

A study of the (static) fission-barrier height B" of nuclei is 
motivated by the importance of this quantity for several physi- 
pal phenomena. Many heavy nuclei decay by spontaneous fis- 
sion, and the size of the fission barrier is a measure of stability 
of the nucleus reflected in the spontaneous fission lifetimes of 
these nuclei lit]. The probability for the formation of a super- 
heavy nucleus in a heavy-ion-fusion reaction is also directly 
connected to the height of its fission barrier [2]. The height B^J 
is a decisive quantity in the competition between neutron evap- 
oration and fission of a compound nucleus in the process of its 
cooling. The large sensitivity of the cross section cr for the syn- 
thesis of the fissioning nuclei on the barrier height B-^J stresses 
a need for accurate calculations of this value. For example, a 
change of B^' by 1 MeV changes the calculated survival prob- 
ability of a synthesized nucleus by about one order of magni- 
tude or even more iQ] ■ The population and survival of hyper- 
deformed states at high spin also depends on the fission barriers, 
see e.g. Refs. flS]- In addition, the r-process of stellar nu- 
cleosynthesis depends (among other quantities such as masses 
and y6-decay rates) on the fission barriers of very neutron-rich 
nuclei ii let]. 

The fission barriers for Actinides and super-heavy nuclei 
show significant differences when calculated by various the- 
oretical approaches (see, for example. Fig. 25 in Ref. Ill 
and Fig. 2 in Ref. f2(]). Many extensive calculations of fis- 
sion barriers have been done in the framework of the macro- 
scopic -i-microscopic (MM) method (see Refs. 
IT2] and references therein). In addition, after the pioneering 
work of Flocard et al. 11311 a growing number of self-consistent 
investigations have been reported in recent years. These calcu- 



lations are based on modern energy density functionals. After 
semiclassical investigations [14] there are now fully^ quantum 



mechanical calculations available based on Skyrme [ 151 1161117 , 
M19,M, Gogny IH tH or RMF 123* 21, i4 iJ function- 
als. The majority of these calculations are restricted to axially 
symmetric shapes. In reality, however, the axially symmetric 
fission barriers sometimes indicate only an upper limit, because 
calculations including triaxiaUty |20l |25l t22l tlSfl have found a 
lowering of the fission barrier by up to several MeV. This low- 
ering strongly depends on the proton and neutron number and 
on the model employed. 

So far all investigations of fission barriers based on covariant 
density functional theory (CDFT) have been performed in the 



RMF+BCS framework ||24|, [15|, [ITD. One of the major goals 



of the current manuscript is to perform the detailed compara- 
tive study of the inner fission barriers in the RMF-hBCS and 
in the relativistic Hartree-Bogoliubov (RHB) frameworks. In 
particular, we try to understand to what extent the properties of 
the first barrier are influenced by the impact of different pairing 
schemes and how the uncertainties in the extrapolation of the 
pairing strength towards super-heavy region affect the fission 
barriers in these nuclei. In order to save computational time we 
do not consider the outer fission barriers in the current investi- 
gation. This restriction has its own merits. The inner barriers 
are generally better measured than the outer ones, and they are 
certainly more important for the r-process, since they determine 
the thresholds. Furthermore, spontaneous fission lifetimes tend 
to be dominated by the inner barrier, even if occasionally an 
outer barrier can have a crucial effect if it is large enough. The 
consideration of only inner fission barriers allows us to restrict 
the calculations to reflection symmetric shapes, because it has 
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been found in several investigations |24, US [H S H that 
odd-multipole deformations (octupole, etc.) do not play a role 
in the inner fission barrier of the Actinides and of superheavy 
nuclei. 

The manuscript is organized as follows. The theoretical 
framework, the selection of pairing schemes and the numeri- 
cal details of calculations are discussed in Sect. |2] In Sect. [3] 
we study the influence of pairing correlations on the shape and 
the height of the inner fission barrier. Experimental data on fis- 
sion barriers in Actinides and superheavy nuclei is compared 
with the results of such calculations in Sect. ID Finally, Sect. |5] 
contains the main conclusions of our work. 

2. Theoretical framework and numerical details 

The calculations discussed in the current manuscript are 
based on two mean field methods to treat pairing correlations 
in nuclei: 



2.1. RMF+BCS calculations 

The RMF-hBCS scheme is rather simple. The RMF- 
equations are solved and at each step of the iteration the BCS 
occupation probabilities are determined. These quantities are 
used in the calculation of densities, energies and new fields for 
the next step of the iteration. Two methods are available for the 
evaluation of the occupation numbers v^, either the "constant 



gap" approximation \ 27] or the solution of the gap equa- 
tions based on a seniority force with the strength parameters 
Gj for neutrons (r - n) and protons (r = p) (denoted in the 
following "constant G"). 

In the first case (constant gap) one starts with fixed gap pa- 
rameters A and uses the BCS expression 



V =r 1- 



Sk- A 



(1) 



with Ek - sjisk - A)^ + A^ for the occupation numbers and for 
the pairing energy 



^pair 



k>0 



UkVk 



(2) 



Ek are the eigenvalues of the Dirac equation and the chemical 
potentials A are determined by the particle numbers. 

In the second case (constant G) one starts with pairing 
strengths parameters G and solves in each step of the iteration 
the gap equation 1128 1 



and determines the gap parameters 



A = G ^ UkVk 

k>0 



(3) 



(4) 



in a self-consistent way. This method is based on the residual 
interaction of the seniority model 1.28.1 and its strength parame- 
ters G. 



In principle both methods are equivalent if the constants G 
are adjusted properly to the corresponding gap parameters A 
(or vise versa). There is, however a difference between the two 
models when we carry out calculations along the fission path by 
constraining the deformation. In the first case (constant gap) the 
gap parameters A are kept fixed as a function of the deforma- 
tion. This means the equivalent strength parameters G change 
with deformation. In the second case (constant G) we use a 
fixed pairing force with constant strength parameters G for all 
deformations and the gap parameters A change with deforma- 
tion. This situation is, of course, closer to reality. 

2.2. Fully self-consistent RHB calculations 

The second method used in this investigation is based 
on a fully self-consistent solution of the relativistic Hartree- 
Bogoliubov (RHB) equations as introduced in Refs. |2^ 30ll : 
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(5) 



where ho is the Dirac hamiltonian of RMF-theory and the pair- 
ing field 



Ai2 = 2vf2V34 



(6) 
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is determined in a self-consistent way by the pairing tensor 

k 

and an effective pairing interaction yP. The pairing energy is 
given by 



= --Tx^K. 



(8) 



In the present investigations we compare two different pairing 
interactions. First we use the Brink-Booker part of the Gogny 
force with finite range 

yPP(l,2)= ^e"<''^'''>'(W,- + BiP^-HiP^-MiP''P^),(9) 

1=1,2 

where P'^ and P^ are the exchange operators for spin and isospin 
and the parameters W/, Bi, //,, and M, (/ = 1,2) of this 
force have been carefully adjusted to the G-matrix calculations 
in nuclear matter and to the pairing properties of finite nuclei all 
over the periodic table (for details of this fit see Refs. ll3lLl32ll '). 
We use here the parameter set DIS 112 ill . 

As second example we use a zero range (J-force 



yPP(l,2) = -yo5(ri -rz). 



(10) 



as it has been used in many non-relativistic HFB calcula- 
tions [18]. This interaction does no depend on density and thus 
it leads to volume pairing. 

In RHB -calculations the pairing gap of BCS-theory is re- 
placed by the pairing field A12 in Eq. (|6]l. In order to have 
a simple measure for the size of pairing correlations in RHB 



2 



theory we therefore introduce in the following calculations the 
"average gap" in the canonical basis 



(A) 



k 



(11) 



where are the eigenvalues of the density matrix = 
^2k^ik ill '^he canonical basis and A^f. are the diagonal ma- 
trix elements of the pairing field A in this basis. For details see 
Ref. Q 



2.3. The pairing window 

Pairing is restricted to the vicinity of the Fermi surface and 
the size of this vicinity is essentially characterized by the gap 
parameter A. This parameter determines the distribution of oc- 
cupations numbers and therefore most of the experimental 
quantities depending on pairing correlations. Of course, in BCS 
or HFB theory based on the seniority model or on zero range 
forces the sums in Eqs. ( I2I3I4I or|7]i show an ultra-violet diver- 
gence and one has to limit the sums in Eqs. ( I2I3I4| | to a pairing 
window Bk < Ecui-oS and in Eq. ^ to Ek < £cut-off- Con- 
cerning this pairing window, there is a difference between the 
methods of constant gap and constant G. In the method of con- 
stant gap, the essential quantity A is fixed and determined in one 
or another way by experiment. The size of the pairing window 
enters only in the calculation of the pairing energy (|2]i which is 
not measurable. For reasonable pairing windows the pairing en- 
ergy is small as compared to the binding energy of the system 
and by this reason one often finds the remark in the literature 
that the results do not depend on the pairing window. However, 
the situation is different for the method of constant G or for zero 
range forces. In this case the size of the effective pairing con- 
stants G or V{) in Eq. (fTOl i have to be adjusted in such a way, 
that the resulting gap parameters A correspond roughly to the 
experimental gaps. As a consequence there is a strong connec- 
tion between the strength of the pairing force and the cut-off 
energy Ecm-oS- It makes no sense to give only the strength pa- 
rameter G (or Vo) if the corresponding value of /icut-oft is not 
known. Essentially one is left with the problem that one has 
one experimental quantity A and two unknowns G (or Vq) and 
^cut-off ■ Usually one chooses a fixed value for iicut-off in a some- 
what arbitrary way and determines the corresponding strength 
parameters in such way that the gap parameter A corresponds 
more or less to the experimental value. 

Bulgac [33] has proposed an elegant way to regularize the di- 
vergence of the gap equation in r-space by removing the diver- 
gent part of the pairing tensor O producing in this way a zero- 
range pairing force, which depends only on one parameter. On 
the other side, as it is clearly seen in momentum space, the two 
parameters G (or Vq) and Ecm-o« are connected to two different 
physical quantities, namely, to the strength and the range of the 
effective pairing force. This means that Ecm-o« should not be 
chosen in an arbitrary way, but it should be determined by the 
range of the force. Of course, it not clear which experimental 
quantity is really sensitive to the range of the effective pairing 
force. To search for such quantities one could do calculations 



with finite range pairing forces of different range adjusting in 
each case the remaining parameters to the experimental gap. In 
this context we will study in Sect. [3]the influence of the pairing 
window of zero-range pairing forces on the size of the fission 
barriers. 

Gogny Jsi , 34 ] has avoided the problem of the pairing win- 
dow by using a finite range force of Gaussian type where the 
range is adjusted to a G-matrix, i.e. via an ab-initio calculation 
directly to the properties of the bare nucleon-nucleon interac- 
tion. The fact that this pairing force is extremely successful in 
reproducing an astonishing number of experimental data both in 

m 



non-relativistic J3 1LI35 
tionals ||37l|38l[l 



and relativistic energy density func- 
shows that this choice of the range is deeply 
connected to the physics of the nuclear many-body problem. 



2.4. Numerical details 



In the current manuscript, we use new versions of the 
RMF-hBCS and RHB codes which were specifically designed 
to describe axially symmetric nuclei with large elongation. 
They allow to use different numbers of integration points in the 
directions along the symmetry axis and in perpendicular direc- 
tion thus allowing a better numerical description of highly elon- 
gated systems. As discussed in Refs. {TL ^ the RMF-i-BCS 
and RHB equations are solved in the basis of an axially de- 
formed harmonic oscillator potential characterized by the de- 
formation parameter ySo and oscillator frequency fiwo = 41A"'^^ 
MeV. The truncation of basis is performed in such a way that 
all states belonging to the shells up to Nf fermionic shells and 
Nb bosonic shells are taken into account. The computational 
time increases considerably with the increase of Nf but it is 
much less dependent on Nb- Thus, a special attention has been 
paid to the selection of Nf of the basis which would allow for 
a systematic study of fission barriers in the nuclei of interest, 
providing at the same time a reasonable numerical accuracy in 
the predictions of the physical observables. 

Extensive tests of numerical convergence have been per- 
formed in the spherical, normal-deformed and superdeformed 
(y?2 ~ 0.7 - 1 .0) minima in the RMF calculations without pairing 
on the example of the nuclei ^^^'U and ^'^"'120 with Z = 120 and 

= 184. Contrary to the previous studies of the convergence 
in the RMF framework which were based on the comparison 
of the Nf and Nf + 2 results, we first defined the "exact" solu- 
tion (extending the calculations up to Nf - 36, Nb = 36) which 
does not change with the increase of Nf, and then found the 
truncation scheme for a basis which provides sufficient numer- 
ical accuracy. It turns out that the binding energies for Nf - 20 
and Nb = 20 are described with an accuracy of 200 keV as 
compared with the exact solution. It was also checked that the 
inner fission barrier is described with an accuracy of ~ 100 keV 
both in the RMFh-BCS and RHB calculations in this truncation 
scheme. Thus all the following calculations are performed with 
Nf^Nb^ 20. 
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Figure 1 : (Color online) Potential energy surfaces in -''"Pu obtained in different 
pairing schemes with the NL3 pai'ameterization of the RMF Lagrangian and 
different values of the scaling factor g. They are normalized to the energy of 
the normal-deformed (ND) minimum. Further details are given in the text. 



3. The role of pairing in defining the inner fission barrier 

3.1. Comparison of different pairing schemes 

The height of the fission barrier depends on the pairing. This 
is a well known fact from the early work of Ref. [41] where it 
was shown for reflection symmetric shapes that barrier heights 
decrease with increasing nuclear pairing. The study of the 
fission barrier in -'*"Pu within the RMF+BCS framework of 
Ref. iH showed that the enhancement of the pairing strength 
of the zero-range pairing force by 20% decreases the fission 
barrier by approximately 2 MeV. Similar results were obtained 
also in non-relativistic HFB calculations based on the Skyrme 
forces in Ref. fis']: the reduction of the pairing strength of zero- 
range pairing forces by approximately 15% leads to an increase 
of the fission barriers by roughly 2 MeV in the Actinides and 
can enhance it by up to 4 MeV in superheavy nuclei thus dou- 
bling the barrier heights. The impact of pairing on fission barri- 
ers has also been studied in the Skyrme-i-HF-i-BCS approach us- 
ing a seniority pairing force and zero-range ^-interactions with 
different forms of density dependence in Ref. | |2o|] . 

In Fig.[T]we investigate as a typical example the inner fission 
baiTier of the nucleus ^'^''Pu first by RMF-hBCS and then by full 
RHB calculations with four different pairing schemes. The ma- 
jor goal of this figure is to show the dependence of the potential 
energy surface (PES) and its profile on the pairing strength as a 
function of the deformation. Therefore we introduce a scaling 
factor g which allows us to change the strength of the pairing 
correlations in an appropriate way. Details will be discussed be- 
low. Note that we have not tried to adjust the pairing strengths 
in all these pairing schemes. As a result, the PES of the dififerent 
panels of Fig. [T] cannot be compared directly. 

Figs. [1^ and b show RMF-hBCS calculations with constant 
gap and with constant G within a pairing window of 

^cut-off ~ 



60 MeV. In the first case shown in Fig. [1^ we use the same 
constant gap parameters for the entire deformation range. 
The size of this pairing gap parameters is determined for the 
ground state in the first minimum by a prescription given in 
Ref.m 

4.8 4.8 
^ MeV, A,^ — 



MeV. 



(12) 



In order to see how the barrier depends in this case on the gap 
parameter we multiply this choice for the gap parameter with a 
scaling factor g 



At ^ gAr 



(13) 



in the range between g - 0.8 and g - 1.2 in a step of 0.1. 
The different potential energy surfaces obtained in this way are 
normalized at the ground state. We observe only small changes 
of the barrier for the different values of g. The height of the 
barrier is the largest for the reduced gap with g = 0.8 and it 
drops only slightly with increasing pairing correlations. 

In Fig. [TJ) we show results of similar calculations, but now 
with constant G. We start in the first minimum with the same 
values for the gap parameters as in Fig. [T^ and determine in a 
self-consistent way for each value of g the strength parameters 



Gr^Ar/J] 



UkVk 



(14) 



In a calculation with constant G this choice leads in the first 
minimum of Fig. [it to the same gap parameters A^ as in Fig.[T^. 
Of course, we also could have scaled the strength parameters 
Gt — > gGj directly. However, because of the strong non- 
linearity of the gap equation this leads to dramatic changes in 
the potential energy surfaces and the two calculations with con- 
stant gap and constant G would have been hard to compare. 
Therefore in Fig. [it we compare the potential energy surfaces 
obtained from constant G-calculations using G^-values derived 
according to Eq. ([T4b after a self-consistent solution for the var- 
ious scaled gap parameters. In this case the calculations for the 
various ^-values in the first minimum are identical in Fig. [T^ 
and Fig.[TJ). However, the remaining part of the energy surfaces 
are very different. In particular the height of the barrier is now 
reduced by nearly a factor 2 and the relative changes for the 
various ^-values are larger than in the case of constant gap. 

This behavior can be understood by the fact that the min- 
ima and maxima (saddle points) along the fission path are in- 
duced by shell effects. The microscopic mechanism of this pro- 
cess is related to the density of the single-particle states in the 
vicinity of the Fermi level which is a function of the deforma- 
tion 144l|45l]. In Fig. [2] we show Nilsson diagrams for protons 
and neutrons in the nucleus ^"^''Pu. Minima in the potential en- 
ergy surfaces are stable configurations. They corresponds to a 
region of low level density, whereas a saddle point occurs in the 
vicinity of level crossings, a region of higher level density. In 
regions of high level density it is easier for the quasi-particles 
to spread around the Fermi surface and therefore the size of the 
pairing correlations depends strongly on the level density. As a 
consequence we find a relatively small pairing gap at the min- 
ima of the PES and large pairing at the saddle points. Of course. 
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Figure 2: Proton and neutron single-particle energies in Pu as a function of the quadrupole deformation /32- They are obtained in the RMF+BCS calculations 
with constant gap approximation and the NL3 parameterization of the RMF Lagrangian. Solid and dashed lines are used for positive and negative parity states, 
respectively. 



this effect is restricted to the cases where we determine the gap 
in a self-consistent way using a pairing interaction, which does 
not change along the fission path. 

Pairing correlations with partial occupations of the different 
levels in the neighborhood of the Fermi surface have the ten- 
dency to wash out shell effects and by this general argument we 
expect a reduction of the barrier height with increasing pairing 
correlations. In fact, we observe this tendency in all four panels 
of Fig. [U However the gap at the Fermi surface is relatively 
small in nuclei (roughly 1 MeV) as compared to the shell ef- 
fects which are of the order of l/ktio {^5.6 MeV in ^"^'^Pu). Thus 
we can understand that in the case of constant gap in Fig. [T^ 
a variation of the scaling factor g in the range of 20% has al- 
most no impact on the PES. This is a new result not available 
in the literature. Note that the constant gap approximation has 
been used in early RMF studies (Refs. |j23l |24t |46t]). In view 
of the results of the present analysis such investigations have to 
be treated with care. On the other hand the calculations with 



Table 1 : Gap parameters for neutron and protons for several deformations along 
the fission path obtained by calculations with constant G„ = 0.0399 MeV and 
Gp = 0.0616 MeV, i.e. for the full black curve in Fig.\^. 
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constant G shown in Fig. [TJ) lead to oscillations of the gap pa- 
rameters along the fission path (see Table[T]i. 

The dramatic change of the fission barrier heights is not 
caused in the first place by the change of pairing correlations in 
total, but rather by the fact that the pairing correlations change 
along the fission path. To elaborate more on this oscillating 
behavior of pairing correlations we show in Fig. [3] the pairing 
energies fipair defined in Eq. ([S) for self-consistent calculations 
with the finite range Gogny-force DIS in the pairing channel. 
We find that these quantities are larger at the saddle point than 
at the normal deformed (ND) minimum. With increasing g, 
the magnitudes of the pairing energies at the saddle point grow 
faster than the ones at ND minimum. As a consequence, ad- 
ditional binding due to pairing (which is related to the pairing 
energies in a highly non-linear way) grows faster at the saddle 
point than at ND minimum with increasing g and as a result, 
the fission barrier becomes lower with increasing g. Note, how- 
ever, that the dramatic changes in the pairing energy cannot be 
seen directly in the change of the barrier height, because they 
are compensated to some extent by the fact that larger pairing 
causes a wider distribution of the occupation probabilities 
around the Fermi surface. 

In Figs. [T]; and [TJl we show potential energy surfaces ob- 
tained from self-consistent RHB calculations based on the pa- 
rameter set NL3 using the Gogny DIS of Eq. (|9]l and the zero 
range 5-force of Eq. ( fTOl i in the pairing channel. For the (5-force 
we use the strength parameter Vo = 300 MeV-fm^ and a pair- 
ing window of /Scut-off - 60 MeV. In both cases we introduce a 
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Figure 3: (Color online) Proton and neutron pairing energies in ^'^''Pu as func- 
tions of /?2-deformation. They are obtained in the RHB calculations with the 
Gogny DIS force for indicated values of g. The arrows show the positions of 
the ND minimum and saddle point. 

scaling parameter of the pairing interaction 

ypp(i,2) gypp(i,2) (15) 

and compare in this way the influence of the strength of the 
pairing force on the fission barrier In both cases we observe a 
considerable reduction of the barrier height with increasing g- 
values. This has the same origin as in the case of the constant G 
calculations in Fig.[T]3. Note, however, that we cannot compare 
directly these results with those of Fig. [TJ), because we do not 
scale the strength parameters Gt of the seniority force by the 
factor g but rather the resulting gap parameters of Fig. [1^. 
This leads to much smaller changes of the parameters Gj- 

Summarizing the results of this section we observe for all 
pairing models in Fig.[T]variations of the PES with pairing. In- 
creasing pairing leads to a reduction of the barrier height. This 
effect is, however, relatively small for calculations with con- 
stant gap. On the other hand for self-consistent calculations the 
oscillations in the level density at the Fermi surface induce pro- 
nounced changes of pairing along the fission path and a consid- 
erable reduction of the barrier heights with increasing strength 
of the effective interaction in the pairing channel. 

3.2. Finite versus zero range pairing 

The basic advantage of the Gogny force is its finite range, 
which automatically guarantees a proper cut-ofF in momentum 
space [47]. Calculations with interactions of finite range re- 
quire a substantial numerical effort and therefore many modern 
energy density functionals like Skyrme functionals are based 
on forces with zero range. In the /^/z-channel, where either only 
the levels up to the Fermi surface are occupied or where the an- 
alytical form of the occupation numbers in Eq. ([TJ leads to 
a fast convergence in momentum space, a gradient expansion 
of finite range interactions leading to zero-range forces of the 
Skyrme type ll48ll are well justified and successful. This is no 
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Figure 4: (Color online) The dependence of the fission barrier height on the 
pairing gap in the ground state shown for the RHB calculations with Gogny 
DIS force and (5-force. The results for (5-force are displayed for 3 different 
values of the cut-off energy i?cui-off- 

longer true in the ;?p-channel, where the analytical form of the 
factors UkVk in the pairing tensor k includes high momenta and 
leads to a ultra-violet divergence. The same is true for the se- 
niority pairing which can be considered as the 7 = part of the 
surface delta interaction |28]. In all these cases one is forced to 
use a cut-ofF energy /icut-oif in order to avoid divergencies. 

The strength parameters of the pairing force Vo (or G) are 
usually deduced for a fixed value of the cut-off energy ficut-off 
from the pairing gaps A extracted from experimental data such 
as the odd-even mass differences at the ground state in the ND 
minimum and it is assumed that the strength parameters do not 
depend on the nuclear deformation. Thus, it is interesting to 
see how the description of the fission barriers differs in the 
RHB calculations with zero range and with finite range pairing 
forces under the condition that the pairing strength parameters 
are defined by the same set of experimental data at the normal- 
deformed minimum. In particular we will investigate whether 
the barrier depends on the choice of the cut-off energy in the 
case of zero range pairing forces. 

In Fig.|4]we study the dependence of the height of the fission 
barrier in the nucleus ^""'Pu on the size of pairing correlations 
in the ground state of the ND minimum. For this purpose we 
carry out RHB calculations with the Gogny force DIS and a 
zero-range 5-force in the pairing channel for various values of 
the scaling parameter g in Eq. ( fTSl l. The resulting barrier height 
is plotted as a function of the resulting average pairing gap ( fTTT i 
in the ground state Ags = j{{A}„ + (A)^). The full (black) curve 
shows the results obtained with the Gogny force DIS and the 
other three lines correspond to calculations with the zero range 
5-force in Eq. ( fTOl l using different cut-off energies /icut-off = 60, 
90, and 120 MeV. This figure clearly indicates that there is a 
strong dependence of the fission barrier height on the treatment 
of pairing correlations. For the Gogny force DIS we find a 
nearly linear decrease of the fission barrier height by roughly 
1.00 MeV per 0.4 MeV change in the gap Ags. For the same 
value of Ags the fission barrier is smaller in the calculations with 
zero range force as compared with the ones based on the DIS 
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Figure 5: (Color online) The dependence of the fission hairier on the cut-oft' 
energy £ciit-off in th£ RHB calculations with (5-force given for 3 different values 
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Figure 6: (Color online) Potential energy surfaces in the Z = 120, W = 172 
nucleus obtained in the RHB calculations with the DD-ME2 parameterization 
of the RMF Lagrangian and different values of the scaling factor g. 



Gogny force. Again we see a decreasing barrier height with 
increasing Agj. However, in the case of zero range force the 
results also depend on the cut-off energy E^ui-os- Only for very 
small pairing correlations Ags ~ 0.5 MeV we have have more or 
less the same barrier heights in the calculations with zero range. 
For larger pairing we observe increasing differences between 
barrier heights calculated with the various values of iScut-oft- 

In particular, for small cut-off energies we find considerably 
lower fission barriers than those obtained with the finite range 
Gogny force. This is shown in detail in Fig. |5] where for the 
same calculations the height of the fission barrier is plotted as 
a function of the cut-off energy /icut-ofF for three values of the 
gap parameter in the ground state (Ags = 0.9, 1.0, and 1.1). 
The dependence of the barrier height on the cut-off energy is 
not eliminated even for these large values of ficut-oif which are 
much larger than the typical ones used in many calculations 
with <5-force (see, for example, Ref. jlm ). 



3.3. Extrapolation to superheavy nuclei 

Systematic experimental spectroscopic data, such as odd- 
even mass differences and the moments of inertia, which allow 
to extract the information on the strength of pairing correlations 
are available only up to the proton number Z » 102 and neutron 
number w 158. With increasing proton and neutron numbers 
such data become scarce and less reliable. Thus, existing pre- 
dictions for superheavy nuclei centered around Z = 120 and 
126 are based on drastic extrapolations involving the changes 
of proton number by more than 20 particles and neutron num- 
bers by 14-26 particles. As with any extrapolation there is a 
considerable degree of uncertainty, and modifications of the 
strength of the pairing interaction by +10% cannot be ruled 
out. Indeed, while providing the average description of pair- 
ing properties modern pairing forces show sometimes apprecia- 
ble local deviations from experiment for physical observables 
which are strongly affected by pairing (see, for example. Fig. 1 
in Ref. 142] and Ref. |49]). 

Fig. |6] shows that modifications of the strength of the pair- 
ing interaction by ±10% have a profound effect on properties 



of superheavy nuclei. For example, a decrease of this strength 
by 10% (the g - 0.9 curve in Fig.|6ll increases the inner fission 
barrier by ^ 1.0 MeV, and thus increases the stability of this 
nucleus against fission considerably. On the contrary, the in- 
crease of the strength of the pairing force by ~10% will modify 
the situation drastically (the ^ = 1.1 curve in Fig.|6l). Indeed, 
we find in this case that the superdeformed minimum becomes 
lower than the spherical minimum. This could indicate that the 
spherical minimum would no longer be the ground state of the 
system since the superdeformed minimum is lower in energy. 
However, one should remember that odd-multipole deforma- 
tions (not included in the cuiTent calculations) are important 
at the outer fission barrier, and thus if included they will most 
likely eliminate this barrier UtIi . As a consequence the second 
minimum will not really stay as a minimum in the PES. In this 
situation, the ground state in the spherical minimum is stabi- 
lized only by the inner fission barrier of the ^ 5 MeV height. 



4. Fission barriers in Actinides and superheavy nuclei 

In Figs. |2] and [8] we compare the results of RHB calculations 
with the parameter sets NL3 |50] and DD-ME2 IH of the RMF 
Lagrangian with available experimental data on the heights of 
the inner fission barrier in Actinides and superheavy nuclei. In 
these calculations the original Gogny force is used (scaling fac- 
tor g = 1.0). It was tested that with this value of g the mo- 
ments of inertia of even-even nuclei, which are very sensitive 
to the strength of pairing correlations, are reasonably well de- 
scribed in cranked RHB calculations with the NL3 force in the 
Nobelium region ifsill and in the lighter Actinide nuclei ifssll . 

Fig.|7]compares the calculated heights E of the inner barrier 
with experiments for Actinide nuclei. While the dependence 
of E on the neutron number is reasonably well described in 
both parameterizations, the dependence on the proton number 
Z is not reproduced. The calculated values of E increase with 
increasing Z while the experimental inner baiTiers almost do 
not depend on Z. The inner fission baiTiers obtained in the cal- 
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Figure 7: (Color online) Experimental and calculated heights of the inner fis- 
sion barriers in Actinide nuclei. Experimental data are taJcen from Table IV in 
Ref. 1 18]. A typical uncertainty in the experimental values, as suggested by 
the differences among various compilations, is of the order of ±0.5 MeV [l&l . 
Results of RHB calculations with the parameter sets NL3 and DD-ME2 of 
the RMF Lagrangian are presented. The dotted lines show the results of the 
RMF+BCS calculations of Ref. flTll which are performed in steps of AN = 4. 



culations with the DD-ME2 parameterization exceed the ones 
obtained with the NL3 parameterization by 1-2 MeV. 

Fig. [8] clearly shows that the fission barriers in superheavy 
nuclei are well described in the RHB calculations with the DD- 
ME2 parameterization. The calculated barrier heights are either 
close to the experimental lower limit or slightly exceed it. Only 
in the case of the (Z = 1 14, = 176) system, the calculated fis- 
sion barrier is slightly below the experimental value. However, 
considering the expected error bars on the height of the fission 
barrier this fact is not important. On the contrary, the results of 
the RHB calculations with the NL3 force considerably (by ~ 2 
MeV) underestimate the height of the fission barriers. 

The RHB results obtained with the NL3 force are very close 
to the ones obtained earlier in the RMF-i-BCS framework with 
the same NL3 force in Ref. [ITJ, see Fig. [T] and Fig. [8] They 
are typically within 0.5 MeV of each other The RMF+BCS 
barriers of Ref. (v^ are either very close to the RHB barriers or 
slightly higher (by several hundreds keV) in the Actinide nuclei. 
On the contrary, they are lower by 0.5-0.8 MeV than the RHB 
barriers in superheavy nuclei. Thus, in going from Actinides 
to superheavy nuclei, the change in the height of the fission 
barriers is more pronounced (by 1 MeV) in the RMF+BCS 
framework as compared with the RHB. 

Based on the RHB results, one can conclude that the calcula- 
tions based on the NL3 parameterization do not leave the room 
for triaxiality and the correlations beyond mean field; this con- 
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Figure 8: (Color online) The same as in Fig.]?] but for superheavy nuclei. Ex- 
perimental data, which specifies the lower limit of the heights of fission barriers, 
are taken from Ref. 1 2] . The assignment of the neutron number has some un- 
certainty, therefore, the same experimental barrier appears for the nuclei with 
the same proton number Z. 



Table 2: Excitation energies E (in MeV) of fission isomers as obtained in the 
RHB calculations with the Gogny DIS force and indicated parameterizations 
of the RMF Lagrangian. Last column shows experimental data taken from 
Ref. lH. 



Nucleus 


E(NL3) 


E(DD-ME2) 


E(exp) 


236u 


1.72 


2.24 


2.75 


238u 


1.70 


1.54 


2.557 


240pu 


2.29 


2.21 


~ 2.8 



elusion is consistent with the one obtained in the RMF+BCS 
framework earlier fr?]. The RHB+NL3 calculations underesti- 
mate the experimental fission barriers in some nuclei already on 
the mean field level. As discussed in the introduction, the triaxi- 
ality can lower the fission barriers by up to 3 MeV or sometimes 
even by larger amount. In addition, the correlations beyond 
mean field can lower fission barrier by as much as 1 MeV UtI] 
(see also Ref. |@]). It is reasonable to expect that with these 
degrees of freedom included the majority of the results with 
the NL3 force will fall below the experiment. On the contrary, 
the experimental fission barriers are overestimated in the RHB 
calculations with DD-ME2 force which indicates that there is 
sufficient room for triaxiality and the correlations beyond mean 
field. It is also interesting to mention that the heights of fis- 
sion barriers as obtained with the density-dependent DD-ME2 
force are close to the typical values obtained in the calculations 
with Skyrme force [17i]. On the contrary, the NL3 parameteriza- 
tion of the RMF Lagrangian, which has no density dependence 
in the isovector channel, systematically predicts lower barriers 
than most Skyrme forces [17]. 

An additional constraint on the selection of the force may 
be provided by the energy of the 0^ fission isomeric state. It 
turns out that for Actinides this energy is less dependent on the 
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variation of the pairing strength than the height of the fission 
barrier (see Fig.[T]l thus providing a more robust probe. Reliable 
values for the energies of the fission isomers in even-even nuclei 
are available only for three nuclei (see Table|2]i. It turns out that 
in average the NL3 and DD-ME2 parameterizations of the RMF 
Lagrangian provide similar results for the energies of the fission 
isomers: these energies are underestimated by at least 0.5 MeV. 

5. Conclusions 

The current investigation of fission barriers within covariant 
density functional theory clearly indicates that when aiming at 
a quantitative understanding of fission properties of heavy and 
superheavy elements, it is important to keep the pairing chan- 
nel under control. It shows that the frequently used constant 
gap approximation provides unphysical results for the height of 
fission barrier. On the contrary, the constant strength approxi- 
mation shows a similar functional dependence of the height and 
the shape of barrier as the RHB approach. Therefore it is im- 
portant to treat the pairing properties in a self-consistent way. 

Seniority forces or zero range forces depend at least on 
two parameters, the strength and the cut-off energy. They are 
usually adjusted to the experimental pairing properties in the 
ground state for a fixed value of the cut-off energy. It is shown 
that this procedure leaves room for uncertainties in the fission 
barriers of up to 1 MeV. 

The extrapolation to superheavy nuclei leads to some uncer- 
tainties in the definition of the pairing strength in the Z = 120 
and 126 regions of superheavy nuclei. It was shown that 10% 
uncertainties in the pairing strength have drastic consequences 
for the structure of these superheavy nuclei. 

The DD-ME2 force is the only presently known parameteri- 
zation of the RMF Lagrangian which provides a good descrip- 
tion of fission barriers in superheavy nuclei in axially symmet- 
ric calculations. Previous and current studies within covari- 
ant density functional theory show that other RMF forces ei- 
ther considerably underestimate the barriers or do not leave the 
room for the triaxiality at the saddle point. 
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